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ABSTRACT 

We study compressible MHD turbulence, which holds key to many astrophysical processes, including star 
formation and cosmic ray propagation. To account for the variations of the magnetic field in the strongly turbu- 
lent fluid we use wavelet decomposition of the turbulent velocity field into Alfven, slow and fast modes, which 
presents an extension of the Cho & Lazarian (2003) decomposition approach based on Fourier transforms. 
The wavelets allow to follow the variations of the local direction of magnetic field and therefore improve the 
quality of the decomposition compared to the Fourier transforms which are done in the mean field reference 
frame. For each resulting component we calculate spectra and two-point statistics such as longitudinal and 
transverse structure functions, as well as, higher order intermittency statistics. In addition, we perform the 
Helmholtz-Hodge decomposition of the velocity field into the incompressible and compressible parts and ana- 
lyze these components. We find that the turbulence intermittency is different for different components and we 
show that the intermittency statistics depend on whether the phenomenon was studied in the global reference 
frame related to the mean magnetic field or it was studied in the frame defined by the local magnetic field. The 
dependencies of the measures we obtained are different for different components of velocity, for instance, we 
show that while the Alfven mode intermittency changes marginally with the Mach number the intermittency of 
the fast mode is substantially affected by the change. 
Subject headings: ISM: structure — MHD — turbulence 
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1. INTRODUCTION 

Astrophysical fluids are magnetized and therefore the as- 
trophysical turbulence is magnetohydrodynamic (MHD) in 
its nature. Compressible MHD turbulence is a ke y ele- 
ment for understan ding star formation (see [M ac Low! 120041 
lElmegreen & Scaloll2004tlMcKee & Ostrikej2007 , and refer- 
ences therein) and velocity fluctuations determine many of its 
properties. For instance, in the modern paradigm of star for- 
mation that turbulent velocity sweep up the matter from large 
expands of the interstellar space to create molecular clouds. 
Thus, it is important to know the statistical properties of the 
velocity field, e.g. its spectrum that reflects how much en- 
ergy is associated with the motions at a particular scale (see 
below). 

A further insight into the properties of turbulence, includ- 
ing its generation, consequences and dissipation calls for the 
use of more sophisticated measures. For example, the pro- 
cesses of magnetic field generation depend on the velocity 
field vorticity associated with the solenoidal motions, while 
the processes of compressing gas are determined by the com- 
pressible motions. In approaching the problem of decom- 
posing velocity field into solenoidal and dilatational part fol- 
lowing the Helmholtz-Hodge decomposition is frequently at- 
tempted (see Federrath et al. 2008, 2009). Another approach 
is the decomposition of the turbulent field in MHD case into 
basic modes, i.e., Alfven, slow and fast waves. While this 
approach is trivial for the case of the strong magnetic fiel d 
with infinitesimal fluctuations (see iDobrowolnv et al.lll980l) . 

(kowaT@ astro .wise . edu 
Jazarian@astro.wisc.edu] 



iCho & Lazarian I (120021 120031) proposed a statistical decom- 
position of modes in the Fourier space. The statistical nature 
of the procedure is clear when one considers its application to 
strongly perturbed magnetic fields. As the Fourier transform 
is defined in the reference frame related to the mean magnetic 
field, while the MHD motions happen in respect to the local 
magnetic field, there is an inevitable contribution of all types 
of motion to the decomposed modes. However, studying the 
cases when the real space decomposition was possible in real 
space, Cho & L azarian I (120031, henceforth CL03) showed that 
the cross-talk between the modes is small for subAlfvenic tur- 
bulence. 

Testing of the results in CL03 and increasing the accuracy 
of the MHD mode decomposition of turbulence is one of the 
key goals of the present study. In doing so in the present paper 
we make use of the wavelet transformations in addition to t he 
Fourier transformations. Wavelets (see Me neveau 111991 alibi) 
present a natural way of describing M HD turbulence. In- 
deed, while in the representation of the iGoldreich & SidhaFI 
(119951 henceforth GS95) model of turbulence the anisotropy 
is frequently described in terms of eddies with parallel and 
perpendicular k± wave vectors, the actual description calls 
for choosing for || and _L in respect to the local magnetic 
fielcfl The latter is really easy to understand, as it is only 
local magnetic field that influences fluid motions at a given 

' Due to this fact that closure relations used for the model justification 
in GS95 are doubtful. The importance of the local system of reference 
was clearly stressed in the works that followed the original GS9 5 study 
rCaz arian & Vishniac 1999; Cho & Vishniac 2000; Maron & Goldreich 
I255lt ICho et aUl200a |2003|: iLithwick & Goldreichl2001l; iCho & Lazarian 
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point. Wavelets allow for a local description of the magne- 
tized turbulent eddies. 

In the paper we decompose the turbulent velocity fields us- 
ing both wavelets and a more traditional Helmholtz-Hodge 
decomposition into solenoidal and compressible parts. We 
feel that the latter decomposition is more justified for the hy- 
drodynamic turbulence than for the MHD turbulence that we 
study here. However, we feel that the use of the Helmholtz- 
Hodge decomposition provides an additional, although lim- 
ited, insight into the properties of compressible motions. 

The three major properties of the velocity field that we 
focus our attention in the paper are turbulence spectra, 
anisotropics and intermittency. These three measures require 
further description that we provide below. 

While turbulence is an extremely complex chaotic non- 
linear phenomeno n, it allows for a remarkably simple statisti- 
cal description (see Biskamp'2003'). If the injections and sinks 
of the energy are correctly identified, we can describe turbu- 
lence for arbitrary Re and Rm. The simplest description of 
the complex spatial variations of any physical variable, X(r), 
is related to the amount of change of X between points sep- 
arated by a chosen displacement 1, averaged over the entire 
volume of interest. Usually the result is given in terms of the 
Fourier transform of this average, with the displacement 1 be- 
ing replaced by the wave number k parallel to 1 and |k| = 1/|1|. 
For example, for isotropic turbulence the kinetic energy spec- 
trum, E(k)dk, characterizes how much energy resides at the 
interval k,k+dk. At some large scale L (i.e., small k), one ex- 
pects to observe features reflecting energy injection. At small 
scales, energy dissipation should be seen. Between these two 
scales we expect to see a self-similar power-law scaling re- 
flecting the process of non-linear energy transfer We shall 
attempt to get the power-law scalings for the components of 
the velocity field. 

The presence of a magnetic field makes MHD tur- 
bulence anisotropic ( Montaomerv & Turner 1981 



Matthae us et"an 1198 31: [ShebaUn et al. 1983; Higdon 1984; 
Goldrei ch & Sidhar ~[l995|)(see lOughton et al. '20031 for 
review). The relative importance of hydrodynamic and 
magnetic forces changes with scale, so the anisotropy of 
MHD turbulence does too. Many astrophysical results, 
e.g. the dynamics of dust, scattering and acceleration of 
energetic particles, thermal conduction, can be obtained if 
the turbulence spect rum and its anisotropy are known (see 
iLazarian et al.l 120091 for review). The knowledge of the 
anisotropy of Alfvenic mode of MHD turbulence deter- 
mines the ext end of magnetic fie ld wandering influencing 
heat transfer (iNarayan & Medvede v 2001; Lazarian 2006) 
and magnetic r econnection Eazarian & Vishniaa il999i ; 
iKowal et alJl2009h . 

We would like to stress that in what follows we discuss 
the properties of strong MHD turbulence. This type of tur- 
bulence is not directly related to the amplitude of the mag- 
netic perturbations, however. The low-amplitude turbulence 
can be strong and isotropically driven turbulence with SB 
at the injection scale exhi bits only a lim ited range of scales 
for which it is weak (see iGaltier et alj 2000), while at suf- 
ficiently small scales it gets strong (see the discussion in 
[Lazarian & Vishniac 1999). 

An anisotropic spectrum alone, say E(k)dk, cannot char- 
acterize MHD turbulence in all its complexity because it in- 
volves only the averaged energy in motions along a particu- 
lar direction. To have a full statistical description, one needs 
to know not only the averaged spectrum of a physical vari- 



able but higher orders as well. The tendency of fluctuations 
to become relatively more violent but increasingly sparse in 
time and space as the scales decreases, so that their influ- 
ence remains appreciable, is called intermittency. The inter- 
mittency increases with the ratio of the size scales of injec- 
tion and dissipation of energy, so the very limited range of 
scales within numerical simulations may fail to reflect the ac- 
tual small scale processes. The turbulence intermittency can 
result in an important intermittent heating of the interstellar 
medium (Falgar one etZll200ll20Ml2007h . 

In this article we investigate the scaling properties of the 
structure functions of velocity and its components for com- 
pressible MHD turbulence with different sonic and Alfvenic 
Mach numbers. In ^we describe the numerical models of 
compressible MHD turbulence. We decompose velocity into 
a set of components including the incompressible and com- 
pressible parts, and MHD waves: Alfven, slow and fast using 
methods described in |3] In ^ we study spectra of velocity 
and its components. In ^we study the anisotropy of dissi- 
pative structures. We show differences in the structures for 
different components. In ^we study the scaling exponents 
and the intermittency of velocity structures. We show their 
dependence on the sonic and Alfven regime of turbulence. In 
we discuss our results and their relation to the previous 
studies. In ^we draw our conclusions. 

2. NUMERICAL SIMULATIONS 

We used an second-o rder-accurate essentially nonoscilla- 
tory (ENO) scheme (see ICho & Lazarian ir2002h to solve the 
ideal isothermal MHD equations in a periodic box. 
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where p is density, v is velocity, B is magnetic field, and a is 
the isothermal speed of sound. We incorporated the fie ld in- 
terpo lated constrained transport (CT) scheme (see, e.g., iTothI 
120001) into the integration of Eq. Q to maintain the V • B = 
constraint numerically. On the right-hand side, the source 
term / represents a random solenoidal large-scale driving 
force. The rms velocity Sv is maintained to be approximately 
unity, so that v can be viewed as the velocity measured in 
units of the rms velocity of the system and B / (AttpY^^ as the 
Alfven velocity in the same units. The time t is in units of the 
large eddy turnover time (^ L/5v) and the length in units of 
L, the scale of the energy injection. The magnetic field con- 
sists of the uniform background field and a fluctuating field: 
B = Bg^i + b. Initially ft = 0. We use units in which the Alfven 
speed va = -Sext/ i^irp)^^^ = 1 and p=\ initially. The values of 
Bext have been chosen to be similar to those observed in the 
ISM turbulence. 

For our calculations, similar to our earlier studies 
( iKowal et all 120071) . we assumed that Bext/(47rp)'''^ 

5B/ (AttpY^^ ~ Sv. In this case, the sound speed is the con- 
trolling parameter, and basically two regimes can exist: su- 
personic and subsonic. Note that within our model, super- 
sonic means low /3, i.e. the magnetic pressure dominates, and 
subsonic means high /3, i.e. the gas pressure dominates. 
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TABLE 1 
List of all analyzed models. 



model 


Bext 


^ini 


Ms 


Ma 


res. 


max time 


max (5p 


max 5V 


BlPl 


1.0 


1.0 


0.68*003 


0.68±003 


512^ 


7.0 


~0.3 


-0.6 


BlP.l 


1.0 


0.1 


2 20^'^-^ 


069±o.03 


512^ 


7.0 


-0.9 


-0.7 


BlP.Ol 


1.0 


0.01 


y.oiO'^ 


0.70±''04 


512^ 


7.0 


-1.9 


-0.7 


B.lPl 


0.1 


1.0 


074±o.06 


7 4±o.6 


512^ 


7.0 


-0.2 


-0.7 


B.lRl 


0.1 


0.1 


2.34±o.08 


7;4±o.3 


512^ 


7.0 


-0.8 


-0.7 


B.lPOl 


0.1 


0.01 






512^ 


7.0 


-1.5 


-0.7 



We present 3D numerical experiments of compressible 
MHD turbulence for a broad range of Mach numbers (0.2 < 
Ms < 7.1 and Ma ^ 0.2 or - 1.8; see TablefB- The model 
name contains two letters: "P" and "B" followed by a num- 
ber The letters B and P mean the external magnetic field and 
the initial gas pressure, respectively, and the numbers desig- 
nate the value of the corresponding quantity. For example, a 
name "B.lP.Ol" points to an experiment with Bgxt = 0.1 and 
P = 0.01. We understand the Mach number to be defined as 
the mean value of the ratio of the absolute value of the local 
velocity v to the local value of the characteristic speed Cj or 
va (for the sonic and Alfvenic Mach number, respectively). 

We drove the turbulence at wave scale k equal to about 2.5 
(2.5 times smaller than the size of the box). This scale de- 
fines the injection scale in our models. We did not set the 
viscosity and diffusion explicitly in our models. The scale at 
which the dissipation starts to act is defined by the numerical 
diffusivity of the scheme. The ENO-type schem es are consid- 
ered to be relatively low diffusion ones (see, e.g. Liu & Osher' 
|T998; Levy, Puppo & Russo 1999). The numerical diffusion 
depends not only on the adopted numerical scheme but also on 
the "smoothness" of the solution, so it changes locally in the 
system. In addition, it is a time-varying quantity. All these 
problems make its estimation very difficult and incompara- 
ble between different applications. However, the dissipation 
scales can be estimated approximately from the velocity spec- 
tra. In the case of our models we estimated the dissipation 
scale ki, at 22 for the resolution 256^. 

3. DECOMPOSITION OF THE VELOCITY FIELD INTO 
COMPONENTS 

3.1. Compressible and incompressible parts 

Using the Hodge generalization of the Helmholtz theorem 
we can split an arbitrary vector field u into three components: 

H=Hp + Hj + H/, (4) 

where each component has specific properties: 

a) Potential component («,,) - it is curl-free component, 
i.e. V X Hp = 0, so it stems from a scalar potential (/>: 

Up = V<t>. (5) 

The scalar potential cj) is not unique. It is defined up to 
a constant. This component describes the compressible 
part of the velocity field. 

b) Solenoidal component (h^) - it is divergence-free com- 
ponent, i.e. V • Mj = 0, so it stems from a vector potential 
$: 

= V X $. (6) 



The vector potential $ also is not unique. It is de- 
fined only up to a gradient field. In the case of velocity 
this component describes the incompressible part of the 
field. 

c) Laplace component (u/) - it is both divergence-free and 
curl-free. Laplace component comes from a scalar po- 
tential which satisfies the Laplace differential equation 
A0 = O. 

Thus the decomposition can be rewritten in the form: 

H = V X $ + V(?!) + H/. (7) 

Applying the divergence operation on both sides and us- 
ing the divergence-free property of the solenoidal and Laplace 
components we obtain 

V-M = VV = A0. (8) 

To find the scalar potential cj) ™d the potential field Up we 
have to solve Poisson equation for V • u. 

To calculate the vector potential $ we apply curl operation 
on both sides of the Eq. dTji. Similarly, using the divergence- 
free property of the potential and Laplace fields we results in 
equation: 

V X M = V X V X $ = A$. (9) 

Here, the calculation of the vector potential $ requires to 
solve triple set of Poisson equations, one equation for each 
component of $. 

The simulations with periodic boundary conditions have 
the advantage, that we can solve the Poisson equation us- 
ing Fourier methods. The Fourier components of the velocity 
field are transformed back to the real space then, and further 
analyzed. 

3.2. Separation into the Alfven, Slow and Fast modes 

Another very important type of decomposition is the sep- 
aration of velocity into the MHD waves: Alfven, slow and 
fast. In this paper, we use an extended mode based on a tech- 
nique described in CL03. The procedure of decomposition 
is performed in the Fourier space by a simple projection of 
the velocity Fourier components u on the direction of the dis- 
placement vector for each mode (see The directions of 
the displacement vectors ^s, £, f, and ^a corresponding to the 
slow mode, fast and Alfven modes, respectively, are defined 
by their unit vectors 

is oc (-1 + a-VD)k\\kii+il + a-VD)k±k± , (10) 
l/cx (-l + a+V^)/t||k||+(l + Q; + VD)/t_Lk_L, (11) 
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real space 



Bo-k plane 



Fourier space 

Fig. 1. — Graphical representation of the mode separation method. We separate the Afven, slow and fast modes by the projection of the velocity Fourier 
component y/^ on the bases ^^i, and ^f, respectively. Figure taken from CL03 
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(12) 



where k|| and kj^ are the parallel and perpendicular to Bgxt 
components of wave vector, respectively, D = (1 + a)^- 
4acos^6', a = a^/V^, 9 is the angle between k and Bgxt, and 
(f is the azimuthal basis in the spherical polar coordinate sys- 
tem. The Fourier components of each mode can be directly 
used to calculate spectra. For other measures, such as struc- 
ture functions, we transform them back to the real space. 

We extend the CL03 technique by introducing an additional 
step before the Fourier separation, in which we decompose 
each component of the velocity field i nto orthogonal wavelets 
using discrete wavelet transform ( see lAntoindl 1 9991 e.g.) 



U{a,wimn) = a ^'^^ ijj 



Xijk - Win 



u(xijk)A''x, (13) 



where Xjjk and wi„i„ are A^-dimensional position and transla- 
tion vectors, respectively, a is the scaling parameter, u{xijk) is 
the velocity vector field in the real space, U{Xijk) is the veloc- 
ity vector field in the wavelet space, and ifj is the ortogonal 
analysing function called wavelet. The sum in the equation 
is taken over all position indices. We use 12-tap Daubechies 
wavelet as an analyzing function and fast discrete version of 
the wavelet transform ( Antoine 1999|), thus, as a result we ob- 
tain a finite number of wavelet coefficients. After the wavelet 
transform of the velocity we calculate the Fourier represen- 
tation of each wavelet coefficient and perform its individual 
separation into the MHD waves in the Fourier space using the 
CL03 method and then update the Fourier coefficients of all 
MHD waves iterating over all wavelets. In this way we obtain 
a Fourier representation of the Alfven, slow and fast waves. 
The final step is the inverse Fourier transform all all wave 
components. 

This additional step allows for important extension of the 
CL03 method, namely, allows for the local definition of the 
mean magnetic field and density used to calculate a and D 
coefficients. Since the individual wavelets are defined lo- 
cally both in the real and Fourier spaces, the averaging of the 
mean field and density is done only within the space of each 
wavelet. 

We can summarize the extended version of the MHD waves 
decomposition into the following steps: 

1. Perform wavelet transform of all velocity components. 

2. Iterate over all wavelet coefficients. 

(a) Calculate Fourier representation of the current 
wavelet. 



(b) Calculate the mean density and magnetic field 
within the space occupied by the wavelet. 

(c) Perform the CL03 separation in the Fourier space. 

(d) Add the contribution from the separated wavelet 
to the Fourier representation of the Alfven, slow 
and fast modes. 

3. Perform the inverse Fourier transform of the Alfven, 
slow and fast wave components. 

4. SPECTRA OF THE VELOCITY COMPONENTS 

In order to obtain one-dimensional (ID) spectra we first cal- 
culate Fourier transform of a quantity and next multiply it by 
its conjunction. We can use this procedure because we used 
a fully periodic domain. The 3D spectra must be averaged 
or integrated over shells k„ <k < k„+i . We used a simple in- 
tegration by summing all squared amplitudes at given shell. 
For each model we have collected data several time steps. We 
used them to increase the size of sample and to measure the 
time departures of the spectra from their mean profiles. We 
should note here, however, that the time averaging and stan- 
dard deviations were calculated in log-log space. Otherwise, 
we could result in taking logarithm of negative number, e.g. 
for ranges of k where the power spectrum has very small val- 
ues (of order 10~* or less). 

In Figure [2] we present spectra for velocity field (top row) 
and its incompressible and compressible parts (middle and 
bottom rows, respectively) for models with different sonic 
Mach numbers. The left column in Figure |2] shows mod- 
els with a strong magnetic field {Ma ^ 0.5-0.7), while the 
right column shows the corresponding models but with a weak 
magnetic field (A4a ^1.7- 1.9). The spectra are averaged 
over several time snapshots taken for a fully developed turbu- 
lence. This allowed us to estimate the spectra time variance 
which is shown as the gray area plotted around the mean pro- 
files. We see that in all cases those variances are small. 

Within the inertial range, which we estimated to be ^ e 
[4,20], the spectra of velocity field slightly change their spec- 
tral indices with the value of Ais- In the case of subAlfvenic 
models the spectral index is close to -5/3. In the case of su- 
perAlfvenic turbulence, the indices do not change with Ms as 
well, and all models show the same, close to -2, value of the 
spectral index at low wave numbers, which is reflection of the 
presence of shocks in the system. Proceeding toward smaller 
scales we observe a bump with spectral index approaching 
-5 /3 which could be explained by the growing importance of 
the mean field at these scales and dominance of Alfven wave. 
The velocity fluctuations have comparable amplitudes for cor- 
responding scales for all models. It means, that the strength of 
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Fig. 2. — Spectra of the velocity field (top row) and its solenoidal and potential parts (middle and bottom rows, respectively) for experiments with different 
sonic Mach numbers for two magnetic regimes: subAlfvenic (left column) and superAlfvenic (right column). The gray area denotes the degree of time variance 
of the spectra. 



fluctuations at particular scale k within the inertial range does 
not depend on A^ ,. 

The total velocity field contains two components: 
solenoidal, which is equivalent to the incompressible part, and 
potential, which contains the compressible part of the field 
and the remaining part which is curl- and divergence-free. In 
Table |2] we show the percentage contribution of each compo- 
nent to the total velocity field. We see, that the compressible 
part constitutes only a fraction of the total field. However, 
the magnitude of this fraction is different for sub- and super- 
sonic models. In the case of subAlfvenic turbulence it is about 
3% in subsonic models and about 7% in supersonic models. 



which confirms a higher efficiency of the compression in the 
presence of supersonic flows. Furthermore, the fraction also 
changes when we compare models with strong and weak mag- 
netic fields. The velocity field, in the presence of a weak mag- 
netic field, contains about 5% of the compressible part in the 
model with M.s ^ 0.7 and even up to 16% in models with 
A^j > 1. The consequence of the presence of a strong mag- 
netic field results in a reduction of the compressible part of 
the velocity field by a factor 2. This indicates a substantial 
role of the magnetic field in the damping of the generation of 
the compressible flows. 
Due to a substantial dominance of the incompressible part 
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TABLE 2 

Percentage amount of the kinetic energy contained within 

EACH velocity COMPONENT. ERRORS CORRESPOND TO A MEASURE 
OF THE TIME VARIATION. 



Ms Ma Vincomp. V^comp. Va Vf 



~0.7 


~0.7 






58±* 


37±3 


4.8±0-' 


~2.2 


--0.7 


93±2 


7±2 


58±5 


33±4 


9±2 


~7.0 


~0.7 


92±2 


7±2 


56±* 


36±4 


8.0±'" 


~0.7 


^lA 


95±2 


5±2 


52*^* 


42±4 




~2.3 






14±2 


47±3 


37±4 


16±2 


~7.1 


~7.1 


84±2 


16±2 


47±4 


33*" 


20±2 



of the velocity field, we expect that its spectra should follow 
the spectra for the total velocity field. In the middle row of 
Figure |2] we present the spectra for the incompressible part 
of the velocity field. We see that the spectra, at least within 
the inertial range, are very similar to those observed for the 
velocity field. Also the spectral indices have very close values 
for models with strong and weak magnetic field. 

In the case of compressible turbulence, more interesting are 
spectra of the potential component (the bottom row in Fig- 
ure |2]). We see that these spectra are different from those for 
the incompressible part. The strength of the compressible part 
confirms the change of contribution of the compressible part. 
It is smaller for subsonic models than for supersonic models, 
which can be seen in the bottom row of Figure |2] However, 
both supersonic models have almost the same spectra with 
the spectral indices about -2 for subAlfvenic turbulence and 
for superalfvenic turbulence. Both supersonic spectra show 
almost exactly the same profiles and amplitudes at all scales, 
even within the dissipation range. This means that the amount 
of compressible part of the velocity field do not change with 
M.S when its value is larger than one. 

The second decomposition, more important for MHD tur- 
bulence, separates the velocity field into three different MHD 
waves: an incompressible Alfven and slow and fast magneto 
acoustic waves, which both are compressible. In Table |2] we 
included the percentage amount of these components in the 
total velocity field. As we see, the most of energy is con- 
tained in the Alfven wave. It is almost 60% in the case of 
subAlfven turbulence, and about 50% for superAlfvenic tur- 
bulence. The slow wave contains approximately 1/3 of the 
total energy. However, for superAlfvenic case, this amount is 
slightly higher. Table|2]suggests, that the slow wave is weaker 
when the turbulence become supersonic. We do not see simi- 
lar behavior for the Alfven wave in the case of models with a 
strong magnetic field. This effect could take place also in the 
superAlfvenic models, but it is weakened by relatively large 
errors. An interesting dependence is observed in the case of 
the fast wave. Although, the fast wave is the weakest among 
all MHD waves, it strongly depends on the regime of turbu- 
lence. Similarly to the compressible part of velocity field, it 
is stronger for models with weak magnetic field. In addition, 
it is much stronger when turbulence are supersonic, but this 
strength seems to be weakly dependent on the sonic Mach 
number. 

In Figure [3] we present spectra for the MHD waves for all 
models from Table [T] In the left column we show spectra 
for models with a strong magnetic field, while in the right 
column spectra for models with a weak magnetic field. Top, 
middle and bottom rows show the Alfven, slow and fast mode 



spectra, respectively. In the case of Alfven wave we see that 
the spectral indices depend on the sonic Mach number weakly. 
All indices lay between -2 and -5/3 with a slight dependence 
on This is similar situation like for spectra of velocity 
field. For superAlfvenic models we observe similar situation, 
however, since the mean field is weaker, the power spectra 
are shifted down to smaller amplitudes. The similarities in 
spectra between the Alfven mode and velocity field are due to 
the fact, that the Alfven mode constitutes the major part of the 
velocity field. 

The slow wave, however, constitutes also a substantial part 
of the velocity field. Looking closer in the spectral indices 
for the slow mode plotted in Figure [3] we see that their values 
do not change significantly with M.s when mean magnetic 
field is strong, although it is difficult to determine one spec- 
tral index do to changing profile of the spectrum with the wave 
number There seems to be stronger dependence of the spec- 
tral indices on the sonic Mach number when the mean field is 
weak, where the inertial region is easier to determine. 

The fast wave, the weakest mode of the velocity field, 
shows two types of spectra depending on the sonic and 
Alfvenic regime of turbulence. In the case of subsonic mod- 
els we see that the fast mode spectrum has index close to -2. 
When the field is weak, the value of index is a bit flatter then 
-2. This indicates a clear dissimilarity of spectra for sub- 
sonic turbulence, for different strength of the magnetic field. 
In supersonic models, however, the spectral indices are going 
closer to the value -2, independently of the sonic Mach num- 
ber. This signifies a growing role of the shock compression in 
supersonic turbulence. 

5. ANISOTROPY 

If we want to study the anisotropy of turbulent structures 
we need to introduce a reference frame. In the case of mag- 
netized turbulence the reference frame is defined in the nat- 
ural way by the local mean magnetic field. The local mag- 
netic field is computed using the procedure of smoothing by 
a three-dimensional Gaussian profile with the width equal to 
the separation length. Since the volume of smoothing grows 
with the separation length /, the direction of local mean mag- 
netic field might change with I at ar bitrary point. This is an 
extension of procedures employed in lCho & Vishniacl (120001) 
and Cho et al. ( 2002). 

To analyze the anisotropy of different components of the 
velocity field we use the total structure function 

5^(0 =(W), (14) 

where (5v is a total increment 

Ml)=|v(x + l)-v(x)|, (15) 

and (...) denotes an ensemble averaging. For each component 
we evaluate the parallel and perpendicular structure functions 
taking respectively the directions of the separation length 1 
parallel and perpendicular to the local mean magnetic field. 
To show the degree of anisotropy we plot the structure func- 
tion in direction perpendicular to the local mean magnetic 
field as a function of the parallel structure function for cor- 
responding separation lengths. In this way we show also, how 
the anisotropy changes with the scale of structures. 

In Figure |4|we present the degree of anisotropy for the ve- 
locity field (top row) and for the solenoidal (incompressible) 
and potential (compressible) parts of the velocity fields (mid- 
dle and bottom row, respectively). The left and right columns 
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Fig. 3. — Spectra of the Alfven, slow and fast mode (top, middle and bottom rows, respectively) of the velocity field for experiments with different sonic Mach 
numbers in two regimes: subAlfvenic (left column) and superAlfvenic (right column). The gray area denotes the degree of departures at single time snapshots 
from the mean profile. 



correspond to the sub- and superAlfvenic turbulence, respec- 
tively. The gray areas under points show the degree of time 
variances of the structure functions used to plot anisotropies. 
In each plot of Figure H] we also plot lines corresponding 

to the isotropic structure, l \\ ^ l±, and theoretical, Z|| ^ 
jGoldreich & Sidhar |[T995 '. GS95). We see that the structures 
of velocity and its incompressible part show the GS95 type 
anisotropy, for both Alfvenic regimes. Another conclusion 
coming from these plots is that the anisotropy is insensitive on 
the value of the sonic Mach number All curves have the same 
shape in the presented plots. The one noticeable difference 
is that in the case of subAlfvenic turbulence the anisotropy 



slightly changes with the values of which are related to 

the scale of structures. Lower values of ^2 correspond to the 
small-scale structures, larger to the large-scale structures. In 
this description we can suspect, that the small scale structures 
are almost isotropic. The degree of anisotropy grows with the 
scale, and for the large scale structures it exceeds somewhat 
the GS95 anisotropy. In the case of a weak external magnetic 
field the degree of anisotropy is uniform over the whole range 
of scales, but also here we notice very good agreement with 
GS95 anisotropy. The compressible component of the veloc- 
ity field behaves differently. For subAfvenic turbulence it is 
almost perfectly isotropic and independent of the sonic Mach 
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number. In the superAlfvenic turbulence, however, we see, 
that structure of subsonic potential field is isotropic on aver- 
age, although there are relatively large time departures from 
the isotropy. On contrary, in the case of supersonic turbulence, 
the structure of the potential field contains a higher amount of 
anisotropic structures and the degree of anisotropy aims to the 
GS95-law. 

In Figure|5]we show the anisotropy of structures for Alfven, 
slow and fast modes of velocity. The separation of the modes 
is performed with the respect of the local mean magnetic field. 
The first comparison of these plots shows, that the degree of 
anisotropy does not depend on the sonic Mach number for 
the Alfven and slow modes. The anisotropy of fast mode, 
similarly to the case of the compressible part of the velocity, 
depends on the regime of sonic motions. If the fluid motions 
are subsonic the structures of the fast waves is more isotropic. 
When the fluid motions are supersonic, more structures be- 
come anisotropic. We should note here, that speaking about 
the anisotropy of structures we understand a mean anisotropy 
of all structures of the analyzed field. Individual structures in 
the turbulence evolve and yield at different deformations re- 
sulting in a constant change of the degree of the anisotropy. In 
this situation, speaking about an anisotropy, we are relating to 
the mean dominant anisotropy of all structures in the system. 

The two remaining modes, the Alfven and slow, show 
very similar anisotropics to those observed in velocity and 
its incompressible part. We see that basically the degree of 
anisotropy of structures in these two modes is closer to GS95. 
In addition, the presence of a strong magnetic field results in 
a much stronger bending of the curves in right plots of Fig- 
ure |5] signifying changes of the anisotropy with the scale of 
the structure. The curves in plots for superAlfvenic models 

are almost straight and independent of the values of . 



and 



^ + l-(l/3//3forC= 1. 



(18) 



6. SCALING EXPONENTS AND INTERMITTENCY 

Intermittency is an essential property of astrophysical flu- 
ids. As intermittency violates self-similarity of motions, it is 
impossible to naively extrapolate the properties of fluids ob- 
tained computationally with a relatively low resolution to the 
actual astrophysical situations. In astrophysics the intermit- 
tency affects turbulent heating, momentum transfer, interac- 
tion with cosmic rays, radio waves and many more essential 
processes. Physical interpretation of intermittency started af- 
ter the work by Koknogorov, but the firs t successful model 
was presented by iShe & Levegue I (119941). T he scaling re- 
lations suggested by iShe & Levegue I (119941) relate ((p) to 
the scaling of the velocity v/ ~ the energy cascade rate 
tj^ ~ Z"', and the co-dimension of the dissipative structures 
C: 



C(p)=-(l-x) + c(l-(l-x/C)''/s 



(16) 



Parameter C is related to the dimensio n of the dissipative 
struct ures D through relation C = 3 — D (iMiiller & Biskampi 
120001) . In hydrodynamical turbulence, according to the Kol- 
mogorov scaling, we have g = 3 and x = 2/3. Vortex fila- 
ments, which are one-dimensional structures, correspond to 
C = 2 (D = 1). In the MHD turbulence we also observe cur- 
rent sheets, which are two-dimensional dissipative structures 
and correspond to C = 1 {D = 2). For these two types of dis- 
sipative structures we obtain two different scaling relations 
(substituting g = 3 and x = 2/3): 



C(p)=^ + 2 



l-(2/3)''/^ 



for C = 2 



(17) 



Relati on (fTTIl is often called the She & Leveque scaling 
(She & Levequel 119941). while relation (TlSl l. the Miiller- 
Biskamp scaling (Mi iller & Biskamp II2000I). There are theo- 
retical arguments against the model (see iNovikov II 1 9941) . but 
so far the She & Leveque scaling is the best for reproduc- 
ing the intermittency of incompressible hydrodynamic turbu- 
lence. 

Structure functions can be calculated with respect to the 
global or local reference frames. By the scaling exponents 
calculated in the global reference frame we understand the 
scaling exponents calculated from the structure functions av- 
eraged over all directions. In the local reference we distin- 
guish between directions parallel and perpendicular to the lo- 
cal mean magnetic field. In this way we define the reference 
frame locally by the direction of the local magnetic field. 

In Figures |6] and |2] we show scaling exponents for the ve- 
locity and all its parts and waves calculated in the global ref- 
erence frame. In the top left plot of Figure |6] we see that for 
the subAlfvenic turbulence the scaling exponents of velocity 
follow the She-Leveque (S-L) scaling with D= 1. Supported 
by the theoretical considerations we can say that most of the 
dissipative structures are one-dimensional. Even though the 
scalings are not perfectly independent of the value of Ais, 
since we see somewhat lower values of ^ for higher p, the dif- 
ferences between these values for models with different sonic 
Mach numbers are within their error bars, thus it is relatively 
difficult to state that the scalings are completely independent 
or only weakly dependent of the values of Ai ,. Looking in 
the corresponding plot for models with a weak magnetic field 
we clearly see that the spread of curves for different sonic 
Mach numbers is much higher than in the previous case. For 
subsonic model the scaling exponents of velocity follow very 
well the theoretical curve defined by the S-L scaling with pa- 
rameter D corresponding to one-dimensional structures. The 
model with A4s ^ 2.3, however, follows perfectly the S-L 
scaling with D = 2 corresponding to the two-dimensional dis- 
sipative structures. Moreover, models with even higher val- 
ues of the sonic Mach number have the scaling exponents for 
p > 3 somewhat below the S-L scaling with D = 2. These 
observations suggest that the scaling exponents of the veloc- 
ity change with the sonic Mach number but only in the case 
of weak magnetic field turbulence. The presence of a strong 
magnetic field significantly reduces these changes and pre- 
serves the generation of the dissipative structures of higher 
than one dimensions. 

After the decomposition of velocity into its incompressible 
and compressible parts we also calculate their scaling expo- 
nents. In the middle and right columns of Figure |6] we show 
the incompressible and compressible parts of the velocity 
field, respectively. The incompressible part it strong. It con- 
stitutes most of the velocity field thus it is not surprising that 
its scaling exponents are very similar to those observed in ve- 
locity. This is true in the case of subAlfvenic models, because 
all curves in the middle left plot in Figure |6] are tightly cov- 
ering the S-L scaling with D= I. The similarity between the 
velocity and its solenoidal part is also confirmed in the case 
of superAlfvenic models but only for subsonic case, when the 
role of shocks is strongly diminished. Two supersonic mod- 
els show exponents following a scaling more closer to the S-L 
one with D=l, yet still with lower values for p >3. The seal- 
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ings of structure of the compressible part, shown in the right 
column of Figure |6] cannot be compared to any of the theoret- 
ical models. Their scalings signify dissipative structures with 
dimensions higher than two, but the theory may be not appli- 
cable here. In order to explain what these scaling exponents 
represent we should describe what could be a physical pic- 
ture of the compressible part of the velocity field. First of all, 
the compressible part is much weaker then the incompressible 
one, thus the structure functions of higher orders amplify rare 
events in the structure, such as individual regions compressed 
by shocks. This is supported by the huge error bars increasing 
with the value of p signifying that these rare events may have 
poor statistics or can change rapidly contributing differently 
at different times and in different models. In such situations 
the scaling exponents for higher values of p are not reliable. 

Another decomposition, which separate velocity into the 
MHD waves, gives an opportunity to calculate their scaling 
relations as well. In Figure|7]we present the scaling exponents 
for the Alfven, slow and fast waves. The Alfven wave is pre- 
sented in plots in the left column. Comparing these plots with 
the corresponding plots for the solenoidal part of the velocity 
reveals obvious conclusion that the Alfveen wave due to its in- 
compressibility should follow exactly the same scaling as the 
incompressible part. This is visible in these two plots clearly. 
Scaling exponents of the solenoidal part and the Alfven mode 
match nicely in plots for subAlfvenic models with a weaker 
dependency of the sonic Mach number and for superAlfvenic 
models where the dependency of A4s is stronger. Two other 
modes, the slow and fast waves presented in the middle and 
right columns of Figure |7] respectively, show similar scal- 
ing relations. For example, the scaling exponents for slow 
and fast waves in the subAlfvenic turbulence follow the S- 
L scaling with D = 2 and depend on the sonic Mach number 
marginally. This signifies that the two-dimensional structures 
dominate in both components. For superAlfvenic turbulence 
the scalings of these two waves show a somewhat different 
situation. Both components, the slow and fast waves, have 
similar values of scaling exponents (, but they change with 
the sonic Mach number. We see that for subsonic turbulence 
most of dissipative structures of the slow and fast waves is 
one-dimensional. Scaling exponents follow very precisely the 
S-L scaling with D= I. In the case of supersonic turbulence, 
however, these scalings suggest the two-dimensional dissipa- 
tive structures again, similarly to the subAlfvenic turbulence. 
This signifies an important role of the magnetic field in the 
generation of the structures of different dimensions in the sub- 
sonic turbulence. 

The above considerations were carried for the scalings rela- 
tions in the global reference frame, i.e., when we do not take 
the direction of the local magnetic field into account. The 
strong magnetic field is dynamically important in subAlfvenic 
turbulence and can greatly influence the generation of the 
structure, thus the statistical methods which include also its 
local direction in the analysis can give a substantial insight 
into the physics of turbulence with the presence of the mag- 
netic field. In Figure|8]we show scaling exponents calculated 
in the local reference frame for velocity and its incompressible 
and compressible parts. Starting from the subAlfvenic turbu- 
lence, which are shown on two upper rows (for parallel and 
perpendicular directions, respectively), we see that the scal- 
ing exponents are different depending on the direction. For 
example, the structures of velocity in the direction perpendic- 
ular to the local field are one-dimensional, as suggested by 
the plot in the second row on the left in Figure [8j This scal- 



ings are very similar to those calculated in the global refer- 
ence frame. However, the same velocity field in the direction 
along the local magnetic field have different structures with 
higher dimensions than two, according to the corresponding 
plots. Both scalings, in parallel and perpendicular directions, 
depend on the sonic Mach number only marginally. Next, the 
solenoidal component shown in the middle column of figure[8] 
have a similar property. Its structures in the perpendicular 
direction are very similar to that observed in the global ref- 
erence frame while in the parallel direction its structures are 
two-dimensional. Again, the potential component, although 
its scaling relation cannot be explained by the current the- 
oretical models, have scaling relations in the perpendicular 
direction very similar to that in the global reference frame. 
All these observations signify that the dominant structures are 
created in the directions perpendicular to the local field, while 
the parallel structures are less significant and usually have 
more dimensions. We see that in these models the role of 
the magnetic field in generation of the structure is very clear. 
What do we expect when the magnetic field is weaker in such 
turbulence. 

Now, we compare the scaling exponents of MHD waves 
obtained in the local reference frame which are presented in 
Figure |9] The Alfven wave which is incompressible and di- 
rected perpendicularly to the local magnetic field is presented 
in the left column of Figure |9] In the subAlfvenic turbulence 
when the magnetic field is strong we expect that the most of 
the structure of the Alfven mode should be generated in the di- 
rection perpendicular to the local magnetic field. This is con- 
firmed by the corresponding plots in Figure |9] The scaling 
exponents in the perpendicular direction are consistent with 
those in the global reference frame (Figure |7ll. The scalings 
suggest the one-dimensional dissipative structures in all sub- 
Alfvenic models independently of the sonic Mach number. In 
the parallel direction the plot suggests the dissipative struc- 
tures with higher direction but at the same time we see very 
large error bars which signify very poor statistics. It means 
that the structures of the Alfven wave created in the direc- 
tion parallel to the local field are marginal. On the contrary, 
the slow wave has direction parallel to the local mean mag- 
netic field, what signifies that its dissipative structure should 
be more dominant in the parallel direction. Indeed, the slow 
wave in the parallel direction has structures consistent with 
those observed in the global mean magnetic field while in the 
perpendicular direction it shows some random, rare events in 
the structure. 

7. DISCUSSION 

7. 1 . Major Accomplishments and Limitations of the Present 

Study 

In the paper we have introduced a new procedure of decom- 
position of MHD turbulence field into Alfven, slow and fast 
modes which uses wavelets. Compared to the decomposition 
procedure based on Fourier transforms described in CL03, the 
wavelet decomposition is more local, thus it follows better the 
local magnetic field direction in respect to which the decom- 
position into modes takes place. As a result, we expect that 
the wavelet decomposition procedure to be more accurate for 
larger amplitudes of turbulence, i.e. larger perturbations of 
magnetic field. 

Our decomposition of the MHD turbulence confirmed 
the results in CL03 in terms of spectra, namely, that the 
Alfvenic and slow modes are anisotropic and consistent with 
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Fig. 8. — Scaling exponents of the velocity and its solenoldal and potential parts (left, middle and right columns, respectively) calculated in the local reference 
frame. Two rows show the scaUngs exponent at the parallel and perpendicular direction the local mean magnetic field for models with a strong magnetic field. 




Fig. 9. — Scaling exponents of the Alfv6n, slow and fast waves of the velocity (left, middle and right columns, respectively) calculated in the local reference 
frame. Two rows show the scalings exponent at the parallel and perpendicular direction the local mean magnetic field for models with a strong magnetic field. 
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the predictions of the GL95 model of incompressible turbu- 
lence, while the fast modes ar e mostly isotropic and form 
an acoustic turbulence cascade dLithwick & G oldreich"2001': 
ICho & Laz arian 2002). As these results were used in stud- 
ies of cosmic ray scattering and acceleration ( Yan & Lazariaji 
200^ 2004, 2008; Cho & Lazarian 20 05; Brunetti & Lazarian 
20071) as well as charged dust acceleration ( Lazarian & Yan 
2002 lYan & Lazarianll2003l: lYan et al.ll200l iYan.,2009.) this 
is an encouraging development. 

At the same time, the intermittency of the different MHD 
modes were shown to be very different. We clearly see the 
dependence of high order statistics of compressible motions 
on the Mach number. We interpret this dependence as the re- 
sult of shock formation, which eventually changes the nature 
of the compressible motion cascade compared to the CL03 
assumptions. 

The limitations of the present study arise from the yet un- 
clear nat ure of the turbulent ca s cade. For instance, it was 
shown in iBeresnvak & LazarianI (120091) that the degree of lo- 
cality of interactions in hydrodynamic and MHD cascade are 
different. Thus even largest available MHD simulations may 
not present the actual inertial range of the cascade, but the 
measured slope may be strongly affected by the extended 
bottle-neck effect of the simulations. In addition, the limited 
range over which Alfvenic turbulence is weak may exhibit a 
rather different scaling of fast modes as a result of the inter- 
actions of the Alfvenic and fast modes ( Chandranll2005h . 

In addition, within the present study we intentionally do 
not consider the scaling of magnetic perturbations. The ve- 
locity and magnetic perturbations for subAlfvenic turbulence 
show some differences, which are rather difficult to study re- 
liably with the available numerical simulations. These dif- 
ferences are not a part of the GS95 picture, but may reflect 
additional yet unclear properties of the MHD cascade (see 
iMiiller & Biskamp 2000). 

In our study we used only the incompressible driving. 
In the presence of the compressible supersonic driving 
dFederrath et alj|2()09l) the scaling looks differe nt, but the ex- 
istence of the inertial range is then questionable. ' Kritsuk et alj 
(|2009) claimed that combining the compressible and incom- 
pressible driving in the Mach number dependent fashion one 
can obtain a better power-law inertial range. This issue re- 
quires further studies. 

The turbulence driving in our study is balanced, in the sense 
that the energy flows in opposite directions are equal. In the 
presence of sources and sinks of turbulent energy, astrophysi- 
cal turbulence is expected to be imbalanced. Our numerical 
studie s of imbalanced turbulence in Beres nyak & Lazarian I 
( 120101) show that the properties of Alfvenic turbulence 
changes substantially in the presence of imbalance. How- 
ever, the degree of sustainable imbalance in compressible tur- 
bulence is still unclear. One expects the density fluctuation 
in turbulent fluid to reflect the incoming waves, altering the 
imbalance. We believe that in high Mach number fluids the 
imbalance is low due to the existence of substantial density 
contrasts. 

7.2. Astrophysical Implications of the Turbulence 
Anisotropy and Intermittency 

Depending on driving astrophysical turbulence may be sub- 
Alfvenic, if the injection velocity Vl is less than Alfven speed 
Va, Alfvenic, if Vl = Va, and superAlfvenic, if Vl > Va- This 
frequently is also described by the Alfven Mach number Ma = 
Vl/Va- Formally, the GS95 model appUes only to incompress- 



ible motions with Vl = Va, or equivalently Ma = 1. Some 
of the astrophysical applications of the model, indeed, use 
the original form of the theory, which substantially limits the 
applications of the theory (see Narayan & Medvedev 2001). 
However, the model can be easily generalized to cover exten- 
sive ranges of superAlfveni c and subAlfve nic turbulence (see 
iLazarian & V ishniac 1999!; lLazarianll2006h . For subAlfvenic 
turbulence with isotropic driving at the scale L an initial weak 
cascade, in which the parallel scale of motions stays the same 
and the spectrum E{k±) ^ is applicable, transfers to the 
regime of strong turbulence at the scale of LM\, for which 
the GS95 critical balance arguments are applicable. For su- 
perAlfvenic turbulence, while up to the scale LM^ the turbu- 
lence is hydrodynamic, it approaches the GS95-type regime 
for smaller scales. Therefore, the relations obtained for MHD 
turbulence that we have studied above can be generalized for 
cases of different intensity of driving. 

While the GS95 model is a model of incompressible tur- 
bulence, our simulation s confirm the numerical findings in 
ICho & Lazarian I (l2002h and CL03 that the scaling of the 
Alfvenic mode in the compressible turbulence is very simi- 
lar to its scaling in the incompressible case. In particular, the 
GS95 anisotropy of MHD turbulence determines the rate of 
magnetic field wandering which is important for many astro- 
physical processes, including the ubiquitous process of mag- 
netic reconnection (Lazarian & Vishniac 1999). Additional 
implications of magnetic field wandering include the diffu- 
sion ofhea^_and_£osmic rays, MHD acceleration of dust etc. 
(see lLazarian et alj2009l for a review). The wavelet approach 
has the potential of increasing accuracy while studying small- 
scale anisotropy in simulations with strongly perturbed mag- 
netic fields. 

FalgaronelFaD (!20S 120061 Hqo I and collabor ators 
(Hily-Blant & Falgarone 2007; Hilv-B lant et~an 120071 and 
refs. therein) attracted the attention of the interstellar commu- 
nity to the potential important implications of intermittency. 
A small and transient volume with high temperatures or vi- 
olent turbulence can have significant effects on the net rates 
of processes within the ISM. For instance, many interstellar 
chemical reactions (e.g., the strongly endothermic formation 
of CH^) might take place within very intensive intermittent 
vortices. The aforementioned authors claimed the existence 
of the observational evidence for such reactions and heating, 
but a more quantitative approach to the pro blem is possible . 
IBeresnva k^ Lazarian (2007) (see also Laz arian et al.ll20()9h 
used the intermittency scaling and calculated the distribution 
of the dissipation rate in the turbulent volumes. In doing so 
they used the fact that lShe & Levequel ([T994) model of inter- 
mittency corresponds generalized log-Poisson distribution of 
the lo cal dissipation rates Dubrulle (1994); She & Waymir^ 
(119951) . The obtained rates of enhancement were not suffi- 
cient to explain the heating required for inducing interstellar 
chemistry (Beresnyak & Lazarian 2007). The same approach 
was used by Pan et al.i (i200 9) who obtained, however, a differ- 
ent result. We believe that one should distinguish shocks from 
vortical motions while calculating the heating induced by in- 
termittency. Our present study show very different scalings 
relevant to these types of motions. 

7.3. Studies of Compressible MHD Turbulence in 
Astrophysical Context 

Numerous studies of compressible MHD turbulence are 
done in the context of star formation (see reviews by 
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iMac Low et al.& KlesseDl2004l:lMcKee & Ostrikeil2007l and 
references therein). Most of these simulations are focused on 
the large-scale appearances of the turbulence, which is de- 
termined by the turbulent driving and do not exhibit any ex- 
tended inertial range of turbulence. 

Search for the universal relations for c ompres sible turbu- 
lence resulted in the rise of interest to the lFleckl (T^983) idea 
of searching universality not for velocity, but for the combina- 
tion of the velocity and density in the form p'/^v. The numer- 
ical study of hydrodynamic compressible turbulence revealed 
that, indeed, the density modified velocity shows the same 
Kolmog orov scaling both fo r low and high Mach number tur- 
bulence dKritsuk et al.ll2007[). Similar effect was c onfirmed in 
our MHD simulations iKowal & LazarianI (l2007h . However, 
the physical justification of this universality is unclear and 
it may result just from the coincidental compensation of the 
change of velocity and density indexes as shocks develop at 
high Mach number turbulence. 

We report steepening of the spectra of compressible mo- 
tions at high Mach numbers. High resolution hydro simu- 
lations (see Kritsuk et al. 2007) show that the velocity spec- 
trum becomes steeper for high Mach number simulations. 
This corresponds to the obs ervational studies of the s uper- 
sonic velocity tu r bulenc e in iPadoan et al.l (120061 l2009h and 
IChepurnov et alj (12006b . These studies are done with the 
Velocity Channel Analysis (VCA) and Velocity Coordinate 
Spectrum (VCS) techniques, which are theory-motivated and 
tested techniques ([Laz arian & Pogosvan 2000, 2004, 200^ 
I2008t IChepurnov etal.ll2008l) . The application of these tech- 
niques should enhance the range of astrophysical turbulent 
velocity fields that can be studied observationall}|3 It is com- 
paring numerics, observations and theory that the progress in 
understanding of turbulence requires. 

8. SUMMARY 

In this article we presented a new technique of decompo- 
sition of turbulent MHD motions into Alfven, slow and fast 
modes. The technique is based on the use of wavelets, which 
provide a more local decomposition compared to the Fourier 
approach in CL02 and CL03. This enables one to have bet- 
ter accuracy of the decomposition of MHD turbulence into 
fundamental modes for higher amplitude of magnetic pertur- 
bations. 

By applying the wavelet decomposition to the results of 
our simulations of compressible MHD turbulence, we investi- 
gated the scaling properties of velocity in compressible MHD 



turbulence for different sonic and Alfvenic Ma Mach num- 
bers. We analyzed spectra, the anisotropy, scaling exponents 
and intermittency of the total velocity and its components cor- 
responding to the Alfven, slow and fast modes. We found that: 

• The amplitude of velocity fluctuations depends on 
only marginally. The lack of significant dependence 
of the velocity fluctuations is also observed for its in- 
compressible part, as well as for the Alfven and slow 
waves. The compressible part of velocity and the fast 
wave show a dependence on TMj, but only for subsonic 
turbulence. In the case of supersonic models, the fluc- 
tuations of the compressible part and fast mode of the 
velocity have comparable amplitudes. 

• The spectral indices depend on Ma in turbulence with 
a strong magnetic field. In the case of turbulence with 
a weak magnetic field only the indices of spectra of the 
fast wave change between sub- and supersonic models. 
For the other components, the spectral indices do not 
change appreciably with the sonic Mach number While 
our conclusions about spectra of fast modes for sub- 
sonic turbulence agree with the CL03 conclusion about 
the acoustic cascade of these modes, we feel that for 
high Mach number we get the spectrum of shocks. The 
anisotropy of Alfvenic turbulence and slow modes is in 
agreement with GS95 theory for both the cases of high 
and low beta plasmas. The velocity fluctuations of the 
fast modes demonstrate isotropy. 

• In the global reference frame, we observe stronger 
changes of the scaling exponents and intermittency for 
velocity and its all components with M.s in the case of 
turbulence with a weak magnetic field. The intermit- 
tency of structures grows with the values of Al ,. How- 
ever, when the external magnetic field is strong, the in- 
termittency for all components depends on the sonic 
Mach number only marginally. In the local reference 
frame, the scaling exponents turbulence depend on the 
direction with respect to the direction of the local mean 
magnetic field. The dependence is stronger for the sub- 
Alfvenic turbulence. 
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